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Abstract 

Coherent signal containing thermal noise is a mixed state of radi- 
ation. There are two distinct classes of such states, a Gaussian state 
obtained by Glauber-Lachs mixing and a non-Gaussian state obtained 
by the canonical probabilistic mixing of thermal state and coherent 
state. Though both these versions are noise-included signal states, the 
effect of noise is less pronounced in the Glauber-Lachs version. Ef- 
fects of these two distinct ways of noise addition is considered in the 
context of atom-field interaction; in particular, temporal evolution of 
population inversion and atom-field entanglement are studied. Quan- 
tum features like the collapse-revivals in the dynamics of population 
inversion and entanglement are diminished by the presence of thermal 
noise. It is shown that the features lost due to the presence of thermal 
noise are restored by the process of photon-addition. 
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1 Introduction 



Thermal radiation is a source of noise in the context of coherent signal. This 
source of noise is unavoidable at finite temperatures, making it necessary to 
understand its effects on experiments. Thermal radiation is characterized 
by the temperature (T) of the source emitting the thermal photons. The 
density operator for a single mode field in thermal statepp is 
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The average number of thermal photons is n, which is related to the tem- 
perature of the source of thermal photons. The photon-number distribution 
(n\pt\n) is a monotonically decreasing function of n, that is, the maximum 
probability is for the vacuum state. On the other hand, a coherent state of 
radiation field is a pure state which is characterized by a complex number 
a. In the number state basis, the coherent state of amplitude apQ is 



\a) = exp 
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Unlike the thermal state, a coherent state has a single peaked photon- 
number distribution, and the peak of the distribution occurs for n ~ |a| 2 . 
Since the thermal state is a mixed state, except at T = OK, any state of 
light that incorporates thermal noise is also a mixed state. Thermal radi- 
ation can be incorporated into the coherent radiation via either mixing or 
Glauber-Lachs (GL) superposition of coherent state and thermal field. The 
GL version corresponds to the unitarily displaced thermal state (DTS)[2j[3], 

p d =b(a)p t D\a) (3) 

where D(a) = exp[cra^ — a* a] is the displacement operator expressed in terms 
of the creation operator a' and the annihilation operator a of the field. These 
states are the quantum versions of classical channels with additive Gaussian 
noise[4]. The DTS can be charaterized as an intermediate state between 
the thermal state and the coherent state, the two limiting cases of the DTS 
corresponding to n = and a = respectively [5]. The other scheme to 
incorporate thermal noise is mixing, which is relevant if the field is obtained 
by probabilistically choosing photons from a coherent source and a thermal 
source. The resultant state is described by 

Pm = (1 - q)pt + q\a){a\, (4) 
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where q is non-negative and less than unity. Such states are relevant in the 
study of noisy quantum channels that transmit either the coherent state \a) 
with probability q or the thermal state pt with probability 1 — q. In the 
subsequent discussion, the state p m is referred as mixed thermal-coherent 
states (MTCS). 

The DTS and the MTCS, albeit constructed out of the coherent states 
and the thermal states, are not equivalent. The DTS is a Gaussian state 
while the MTCS, which is a mixture of Gaussian states, namely, the coherent 
state and the thermal state, is not a Gaussian state unless q = 0, If 
the mean photon-number n is zero, the DTS becomes the coherent state 
a pure state with no thermal noise; in the same limit, the MTCS is a 
mixture of | a) and the vacuum state 1 0) . It is further required to have q = 1 
to ensure that there is no thermal noise in the MTCS. If the amplitude a 
vanishes, the DTS becomes the thermal state pt while the MTCS is a mixture 
of thermal state and vacuum state. In FigJU photon-number distribution is 
shown for DTS (dashed curve) and MTCS (continuous curve) corresponding 
to a = \/l0 and q = 0.5. The photon number distribution of DTS exhibits 
a single peak, whereas MTCS has one at n = and the location of the other 
peak depends on the values of a and n. 

The paper is organized as follows: in Section II, the DTS is shown to 
be a mixture of photon-added coherent states. In Section III, interaction 
of a two-level atom with the radiation field is studied. The objective of 
the study is to highlight the differences in the evolution of entanglement 
and population inversion when the atom interacts with the DTS and the 
MTCS. In Section IV, effects of photon-addition on atom- field dynamics are 
discussed, followed by a summary of the results. 



2 Glauber-Lachs states as mixture of photon-added 
coherent states 

A photon-added coherent state (PACS) of order m, is defined as[7j 

\a,m) = —= ■ (5) 

^/m\L m (-\a\ z ) 

Here L m (x) is the mth order Laguerre polynomial[8j and m is a nonnegative 
integer. The first index a refers to the complex amplitude and the second 
index m is the order of the PACS. It is important to note that PACS of 
different orders are not orthogonal to each other. These states have been 
studied, both theoretically and experimentally, in a variety of contexts, such 
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Figure 1: Photon-number distribution P(n) for DTS (dashed) and MTCS 
(continuous). The amplitude a is VTO and q = 0.5. 



as nonclassicality and quantum-classical correspondence [9]. In this section 
it is shown that DTS is expressible as a mixture of all orders of PACS, all 
of same amplitude. 

The definition of DTS given in Eq. [3] yields the following expression 
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which is a mixture of orthogonal states, namely, the displaced number states 
|n, a) = D(a)\n). Another way of expressing the displaced thermal states is 
to note that 
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Defining n = l/(exp 7 —1), the above expression for can be simplified to 

p d = ee - e HV aa V 7at V aa . (8) 



The parameter e is given by 1/(1 + n). Expanding e 7 ° a in terms of the 
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Fock state projectors \n)(n\, the density operator of the DTS becomes, 

p d = ee - e H 2 ^- r e eaat a tn |0)(0|a n e M * a . (9) 

n=0 

Since exp(aa^) [respectively, exp(— a*a)] and a) [respectively, a] commute, 
the order of multiplication can be reversed. It is easy to recognize that 
exp(— eaa^)\0) = exp(e 2 |a| 2 /2)|ea), a coherent state of amplitude ea, apart 
from the multiplication factor exp(e 2 |a| 2 /2). Further, a* n \eat) are PACS, 
apart from the normalization factor \Jn\L n {— e|a| 2 ). With these identifica- 
tions, the previous equation is recast as 

n 

L n {— \a\ 2 ) |d, n)(d, n\ , (10) 

where a = a/(l + n). The expression reveals that DTS is a mixture of all 
orders of PACS of amplitude a/(l + n). It is interesting to contrast the two 
expressions for the DTS, given in Eqns. [6l and [TOl respectively: the former 
relation expresses the DTS as a mixture of orthogonal states while the later 
expresses it as a mixture of nonorthogonal states. 
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3 Population inversion and entanglement dynam- 
ics 

In this section, dynamics of a two-level atom interacting with a single mode 
field is considered. The initial state of the field is either DTS or MTCS. 
The presence of thermal noise in the initial state affects the dynamics of 
interaction. Evolution of atom-field state is studied to understand the ef- 
fects of different ways of incorporating the thermal noise. The atom-field 
interaction is modelled by the Jaynes-Cummings hamiltonian|ll| 

Hjc = fiuja)a + ~^ G z + h\(a)a- + aa + ), (11) 

where v is the atomic transition frequency, u) is the field frequency and A 
is the coupling constant for the atom-field interaction. Using a suitable 
interaction picture and rotating wave approximation [12] . the hamiltonian is 

H[ = h — a z + ihX (a)o- — acr+) , (12) 
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where A = uj — v is the detuning parameter. The operators relevant to 
the two- level atom are the Pauli operator a z and, the raising and lowering 
operators cr_ and a + respectively. 

The time-evolved density operator p(t) is p(t) = U (t) p(0)W (t) , where 
p(0) is the initial density operator of the atom-field system. The unitary 
evolution operator U(t) is exp —iHjt/h . In this report, the initial state of 
the atom- field system is taken to be |e)(e| ® pf, a factorizable state which 
has no entanglement. The operator pf stands for the field state which is 
taken to be either DTS or MTCS. The states |e, n) and \g, n + 1) transform 
under the action of Hj as follows: 



Hj\e,n) = H— |e,n) + ih\\/n +T|<7, n + 1), 

A _ 

-ff/lfl 1 , n + 1) = ft— \g, n + 1) — ih\y/n\e, n). 



(13) 
(14) 



The results indicate that the span of the states |e,n) and \g,n + 1), is 
invariant under the action of the unitary evolution operator exp(—itHi/h). 
If the interaction is resonant, (A = 0), the time-evolved density operator is 
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The superscript r implies matrix transpose. The coefficients c nm are the ma- 
trix elements {n\p\m) of the initial density operator of the field; the relevant 
expressions for the DTS and the MTCS are 
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respectively. It is also clear that c mn = c* m , the complex conjugate of c nm . 

For the purpose of comparing the interaction dynamics of the atom with 
the two classes of the field, namely, the DTS and the MTCS, it is essential 
that a parameter is chosen to compare the two fields. While the displacement 
a and mean number of thermal photons n are common to both the states, 
the parameter q in the MTCS is not fixed. A meaningful way is to choose the 
value q so that the overlap with the coherent state |a) is the same for both 
the DTS and the MTCS, i.e., (a\p m \a) = (a\pt\a). With this restriction, the 
value of q is fixed by the parameters a and n. The respective contributions 
of coherent state |a) to the DTS and MCTS are 

1 — q Ittl 2 
(a\p m \a) = q + — — exp[- — — , (19) 

1 + 71 1 + 71 

(a\p d \a) = (20) 

1 + 71 

The contribution from the coherent state \a) to the states decreases as n 
increases. Let q that special value of q which ensures equal overlap of the 
DTS and MTCS with the coherent state \a). Then, the last two expressions 
yield 

__ l-exp[-|q|7(l+n)] 
n + 1 — exp[— |a| 2 /(l + n)] 

In physical terms, this is the probability with which coherent state photon 
is chosen in the mixed state defined in EqU] to ensure equal overlap with 
| a). If n is zero (no thermal photons), q becomes unity. Hence, noise-free 
limit of the MTCS is the pure coherent state \a), as in the case of DTS. If 
a = 0, then q is 0. In this limit, the states pd and p m become the thermal 
state pt, a mixed state. Thus the condition of equal overlap implies that 
the limiting cases of the DTS and MTCS are the same. This implies that 
the DTS and the MTCS are two classes of states that interpolate between 
a pure state (|a))and a mixed state (thermal state). 

During evolution, the state of the atom and the field is entangled, a fea- 
ture that originates in the bipartite nature of the system. The entanglement 
between the atom and the field is quantified by the logarithmic negativity 
N, defined as the absolute sum of the negative eigenvalues of the partially 
transposed density operator /3 pr |13j. If are the eigenvalues of p PT (t), 
then N(t) = J2k[\^k\ — Afc] /2. In the case of atom-field evolution under 
Jaynes-Cummings interaction, the time-averaged entanglement depends on 
the initial mixedness of the bipartite state|14j. The mixedness of the DTS 
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1 -Tr[p 2 d } = ^—. (22) 
lHdS l + 2n y J 

This expression is independent of a and hence the mixedness of the DTS is 
same as that of the thermal state. This is a consequence of the fact that the 
two states are related by a unitary transformation. The mixedness of the 
MTCS, given by 

-i m r»2i i 2 (!-g) 2 2g(l -q) \a\ 2 
l-Tr[ Pm ]=l- q - T -^- n __exp[-— ], (23) 

depends on n and a. 

As far as the atomic dynamics is concerned, the relative probabilities to 
be in the excited and states are important. In this context, A quantitative 
population inversion W(t), defined as the difference of the probabilities for 
the atom to be in the excited state and ground state, is of interest. Since the 
atomic Hamiltonian is proportional to \e) {e\ — \g) (g\ , the population inversion 
is essentially the energy of the atom, apart from an overall scaling factor Tiv. 
If the state of the field is a coherent state of suitably large amplitude, the 
time-evolution of Wit) is characterized by the well known "collapse-revival" 
structures [El EH E] . There are associated features present in the evolution 
of entanglement between the atom and the fieldpQ. The effectiveness of the 
thermal radiation in diminishing these features of atom-field is a qualitative 
indicator of the effects of thermal noise [18] . 

If n and \a\ are less than unity, both the DTS and the MTCS are ap- 
proxiamtely the vacuum state of the field. The interaction dynamics of a 
two-level atom in its excited state with the vacuum state leads to a periodic 
exchange of energy between the field and the atom. Therefore, nontrivial 
dynamics can be expected if |q| differs significantly from unity. As a rep- 
resentative case, the numerical values considered for the study in this work 
are: \a\ 2 = 10 and n is taken to be a fraction of \a\ 2 . 

In Fig. [21 the temporal variation of population inversion is shown for 
three different values of n. For easy comparison, the corresponding profiles 
for the DTS and MTCS are shown in adjacent columns. In the noise-free 
limit, that is, n = 0, the DTS and the MTCS are identical and there is no 
difference in the dynamics of interaction with the two- level atom as seen from 
the first row of figures. However, as n is increased to 0.1, equal to one percent 
of the coherent state photon number, the dynamics in the DTS case is very 
different from that of the MTCS. This is seen by comparing the two figures in 
second row. The profiles clearly show the marked differences in the dynamics 
induced by the two different initial states. As the fraction of thermal photons 
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Figure 2: Time evolution of population inversion W(t). In all the subplots, 
the abscissa corresponds to the scaled time Xt. Left column corresponds to 
evolution with DTS as the initial state of the field and the right column 
corresponds to MTCS as the initial state. For the both the states a = \/l0- 
Different rows correspond to different values of mean number of thermal 
photons: n = (first row), O.l(second) and 1.0 (third). 

is increased to ten percent (n = 1), the temporal profile of the population 
inversion in the case of MTCS exhibits higher degree of deviation from the 
n = case, notably, the collapse-revival structure is totally absent. In 
contrast, in the case of DTS the collapse-revival structure is still intact 
even if h = 1. This indicates that the thermal noise in the MTCS is more 
effective than that in the Glauber-Lachs version. The temporal evolution 
of negativity N(t) is shown in Fig. [3j To compute the eigenvalues of the 
partially transposed bipartite density matrix to high accuracy, the size of 
the density matrix has been chosen to be 300X300. Since the initial photon- 
distribution is appreciable for n ~ 10 and falls rapidly for larger values of 
n, the chosen size is large enough to ensure that the computed eigenvalues 
are not spurious. The first, second and third rows of figures in Fig. |3] 
correspond to n = 0.1, 1 and n = 1 respectively. As in the case of population 
inversion, the dynamics of N(t) shows significant changes in the case of 
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MTCS as the value of n increases, leading to complete loss of structures 
in the temporal profile. The Glauber-Lachs mixing too leads to changes in 
the evolution of N(t) as noise level increases, though to a lesser extent than 
in the MTCS case. Thus, atom-field entanglement is more robust against 
thermal noise if the thermal noise is included through Glauber-Lachs version 
rather than mixing. It may be noted that if \a\ is large, then the condition 
of equal overlap with coherent state also implies that the states are of equal 
mixedness. 



4 Noise reduction by photon-addition 

In the previous section the effect of thermal noise on the dynamics of atom- 
field interaction was discussed. In this section, effects of photon-addition are 
presented. One of the consequences of photon-addition is that the resultant 
state has no overlap with the vacuum state 1 0) . In the thermal state of the 
radiation field, the vacuum state probability is not less than that of other 
number states. Since the photon- added states do not have any contribution 
from the vacuum state, it can be expected that the effects of noise are less 
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pronounced. 

The photon-added versions of the DTS and MTCS are 



Pm = =r-Ff it, 24 

pd = iWp^Y (25) 

wherein the photon-added states are indicated with a tilde. 

Photon-addition introduces nonclassical features such as sub-Poissonion 
statistics, squeezing, etc, apart from increasing the mean number of photons [19 
\20\ l2l] . Additionally, effect of photon-addition shows up in the atom-field 
interaction dynamics too. In particular, the features lost due to the presence 
of thermal photons are restored as a consequence of photon-addition. The 
overlap between the coherent state \a) and the photon-added versions of the 
DTS and MTCS are related to the corresponding overlaps of the DTS and 
MTCS via, 

\a\ 2 

(a\p m \a) = , =rj— i2(a|p m | a), (26) 

1 + qn + (1 — q)\a\ z 

\a\ 2 

(a\p d \a) = 1 + - +N (a\p d \a). (27) 

If |a| 2 = 10 and n = 1, the probability q is nearly 0.5. It is seen that the over- 
lap of the photon-added MTCS p m with the coherent state is significantly 
larger compared to (a\p m \a). So, it is natural to expect the dynamics of 
atom-field interaction to be similar to that in the case of field being in a co- 
herent state. It is interesting to note that in the case of pt, photon-addition 
leads to a decrease in the overlap with coherent state. 

The evolution of population inversion in the case of photon-added states 
shown in FigUJ The effects of photon-addition on the dynamics of W(t) is 
clear on comparing the corresponding profiles in Fig. 0] and Fig. [2j Noise 
does not modify the profiles much if the initial state is the DTS or its photon- 
added version. However, comparing the right columns of the two figures, it 
is clear that for the MTCS, the effect of photon-addition is very significant. 
In particular, with n = 1, the evolution of W(t) in the photon-added version 
is very similar to the zero noise, that is, n = case. 

In Fig. [5l evolution profiles of negativity if the initial state is one of 
the photon-added states are shown. To assess the effect of photon-addition, 
the corresponding profiles in Figs [3] and [5] are compared. Specifically, the 
right columns of the two figures show the recourse to noise-free evolution 
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Figure 4: Time evolution of population inversion in the interaction between 
atom and photon-added states. The abscissa corresponds to the scaled time 
Xt. Left column corresponds to evolution with the photon-added DTS p^ as 
the initial state of the field and the right column corresponds to p m as the 
initial state. For the both the states a = \/l0. Different rows correspond 
to different values of mean number of thermal photons: n = (first row), 
0.1 (second) and 1.0 (third). 



due to photon-addition to the MTCS. As in the case of W(t), the effect 
of photon-addition is very mild in the case of DTS. Though the process of 
photon-addition is applied to both the coherent component and the thermal 
component of the MTCS, the net effect is to produce effects similar to that of 
coherent radiation as far as the atom- field dynamics is concerned. In short, 
the noise effect is suppressed to a large extent resulting in the recovery of 
the collapse-revival structures present in the noise-free evolution. 

The reason for the recovery of noise-free evolution on photon-addition 
is due to the the fact that the vacuum state probability is vanishes. The 
photon-added states p m and pd have no overlap with the vacuum state |0), 
that is, (0|p|0) = 0, where p refers to the photon-added versions of the DTS 
and MTCS. When the mean number of thermal photons is about unity, the 
major contribution to the thermal state is from the vacuum. In the photon- 
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Figure 5: Time evolution of negativity N(t) due to the interaction between 
atom and photon-added states. Other descriptions are same as given in Fig. 
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added versions, the contribution from the vacuum is absent, resulting in 
lesser overlap with the thermal state. This, in turn, means that the overlap 
with coherent states is more for the photon-added versions of the MTCS 
and DTS. 

5 Summary 

Thermal states and coherent states are Gaussian states. While displacing 
the thermal state is a way of incorporating thermal noise, and the resultant 
state is still Gaussian. Mixing of thermal and coherent states is another way 
of including thermal noise, however, the resultant state is non-Gaussian. 
The displaced thermal state is a mixture of photon-added coherent states of 
all orders. Though the two classes of states are closely related to thermal 
and coherent states, it is only for a special choice of mixing probability in 
the mixture of thermal and coherent states, these two classes of states have 
equal overlap with coherent states. In the case of displaced thermal state, 
the dynamics of population inversion and atom-field entanglement are not 
very susceptible to thermal photons if the coherent state amplitude is large 
(a ~ \/l0) and the thermal photons contribution is as high as one tenth 
of the total number of photons. In the case of MTCS, though it has the 
same overlap with the coherent state \a) as with the case of DTS, increasing 
the thermal photon contribution renders the dynamics totally devoid of any 
collapse-revival structure in the evolution of population inversion and atom- 
field entanglement. Photon-added versions of these two classes of states 
exhibit, in the context of atom-field interaction, features very similar to 
the case of dynamics in the absence of thermal noise. In essence, photon- 
addition suppresses the effects of thermal noise, thereby amplifying the effect 
of coherent states on the dynamics. 
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